International Journal of Engineering Inventions 

e-ISSN: 2278-7461, p-ISSN: 2319-6491 
Volume 2, Issue 2 (January 2013) PP: 01-09 



Experimental and Theoretical Study of Heat Transfer by 
Natural Convection of a Heat Sink Used for Cooling of 

Electronic Chip 

'Sunil Hireholi, 2 K.S. Shashishekhar, 3 George. S. Milton 

123 Dept. of Mechanical Engineering 
12 S. I. T, 3 S. S. I. T 
123 Tumkur, Karnataka, India 



Abstract:- In the present work, heat transfer analysis of a commercially available heat sink is carried out. The 
heat sink is used for cooling of electronic chip STK4141II and removes the heat generated by the chip through 
natural convection. The work includes experimental investigations and theoretical modelling. 
In experimental investigations electrical heating is used to supply heat to the heat sink and the temperature of 
the heat sink is measured using RTD thermocouples attached to the heat sink. Theoretical temperature of heat 
sink is predicted employing 2-D modelling of the heat transfer process based on fundamental heat transfer 
principles. Theoretically predicted temperatures of the heat sink are compared with the measured temperatures. 
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I. INTRODUCTION 

Advances in the field of electronics have resulted in significant increase in density integration, clock 
rates and emerging trend of miniaturization of modern electronics. This resulted in dissipation of high heat flux 
at chip level. In order to satisfy the junction temperature requirements in terms of performance and reliability, 
improvements in cooling technologies are required. As a result thermal management is becoming important and 
increasingly critical to the electronics industry. The task of maintaining acceptable junction temperature by 
dissipating the heat from integrated circuit chips is significant challenge to thermal engineers. The electronic 
cooling is viewed in three levels, which are non-separable. First, the maintenance of chip temperature at the 
relatively low level despite the high local heat density. Second, the heat flux must be handled at system or 
module level. Finally, the thermal management of the computer machine room, office space or 
telecommunication enclosure [1]. 

In electronic systems, a heat sink is a passive component that cools a device by dissipating heat into the 
surrounding . Heat sinks are used to cool electronic components such as high-power semiconductor devices, and 
optoelectronic devices such as higher-power lasers and light emitting diodes (LEDs). In the present work we 
carry out analysis of heat transfer by a commercially available heat sink used for chip level cooling. 

II. LITERATURE SURVEY 

The cooling aspects have been studied by many investigators. Among them Mobedi and Sunden [2] 
investigated the steady state conjugate conduction-convection on vertically placed fin arrays with small heat 
source inside. Thermal performance of a free standing fin Structure of copper heat sink has been reviewed by 
Jan Bijanpourian [3]. In his work he had used CFD to evaluate and compare heat dissipation capabilities of 
aluminium and copper heat sinks. Experimental Investigation of Pin Fin Heat Sink Effectiveness for forced 
convection is reviewed by Massimiliano Rizzi et al [4]. One of the earliest studies about natural convection heat 
transfer from fin arrays was conducted by Starner and McManus [5]. Four different fin array configurations with 
three base types were investigated and heat transfer coefficients were calculated. Leung et al. [6] performed an 
experimental study on heat transfer from vertically placed fin arrays produced from an aluminium alloy . It was 
found that for different configurations the maximum heat transfer rate from the fin arrays was obtained at the fin 
spacing value of 10 mm. Harahap et al. [7] performed experiments on miniaturized vertical rectangular fin 
arrays in order to investigate the effect of miniaturizing on steady state rate of natural heat transfer. A numerical 
analysis on natural convection heat transfer from horizontally placed rectangular shrouded fin arrays were 
performed by Yalcin et al. [8]. Commercially available CFD package PHOENICS was used to solve the three 
dimensional elliptic governing equations. Yiincii and Mobedi [9] investigated the three dimensional steady state 
natural convection from horizontally placed longitudinally short rectangular fin arrays numerically. A finite 
difference code in Cartesian coordinate system based on vorticity-vector potential approach was used to solve 
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the problem. Kundu and Das [10] performed an analytical study to investigate performance and optimum 
design analysis of four fin array types viz, longitudinal rectangular fin array, annular rectangular fin array, 
longitudinal trapezoidal fin an-ay and annular trapezoidal fin array under convective cooling conditions are 
investigated. Another experimental study regarding the natural convection heat transfer from both vertically and 
horizontally oriented fin arrays was done by Leung and Probert [11]. The effects of fin spacing and fin height 
were investigated for a limited number of fin array configurations. Vollaro et al. [12] analysed natural 
convection from rectangular and vertical finned plates numerically in order to optimize the fin configuration. 
The maximum heat transfer rate from fin array was investigated for the optimum fins spacing as a function of 
parameters such as dimensions, thermal conductivity, fins absorption coefficient and fluid thermo -physical 
properties. 

III. EXPERIMENTAL STUDIES 

The front, top and side views of the heat sink analysed are shown in Fig. 1. It consists of a vertical base 
and seven horizontal fins as per the dimensions shown in the figure in millimetres. In the experimental studies 
electrical heating is employed to supply heat to the heat sink. The electrical heater has same nominal dimensions 
as the electronic chip and occupies the same physical location on the heat sink as the chip. The thermocouples 
used for experiment were calibrated. 

A. Experiment to determine heat transfer coefficient with heating coil 

The method used for determining the heat transfer coefficient is to give a known heat input generated 
by electrical heating and measure the temperature attained by the heat sink. The temperature is measured using 
the RTD-thermocouples at three locations on the model and average of all the three values, T b T 2 , T 3 , is taken 
for the calculation of overall heat transfer coefficient. It is made sure that steady state conditions are reached 
before the temperature are taken and then the corresponding heat transfer coefficient is computed. The 
temperatures are taken at regular intervals till steady state is reached. The values of the temperature after it has 
reached steady state are used in calculating the corresponding heat transfer co-efficient for the given heat input. 
This is done for varying values of heat inputs 
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Fig.l Front, top and side views of the heat sink 



Fig. 2 shows a photograph of the Experimental set-up used for measurement of temperature. The heat 
sink is kept at a height from the surface of the table to enable free flow of air from the bottom of the sink. 
Alternating current from mains is supplied to the heating coil attached to the heat sink through a dimmerstat 
(auto-transformer) as seen in the photograph given in Fig. 2. Voltage across the heating coil is measured using a 
voltmeter and the current is measured using an ammeter, as described above. Heat dissipated and the associated 
heat transfer co-efficient are calculated using the following basic equations [13], [14]. 



cinput 



'avg 

h=- 



= VI 

T1+T2+T3 



Qinput 



(1) 



A(Tavg -Tamb ) 

Where, 

Qinput = Heat input by the heating 
coil (Watts) 



(2) 
(3) 
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V = Voltmeter reading (Volts), 
I = Ammeter reading (Amperes), 
T 1 ,T 2 ,T 3 =Temperature readings of 

Thermocouples (°C) 
T avg = Average temperature of heat sink (°C), 
h = Heat transfer coefficient (W/m 2 K) , 
A = Effective area of the heat sink available for 

heat transfer, 
T am b= Ambient temperature (°C) 

Results are given in Table 1. 

heat sink 




Fig.2 Experimental setup 

Table 1 shows time in seconds and the corresponding temperature readings of the three thermocouples 
Ti, T 2 , T 3 for different values of the heat inputs along with their voltage and current values. The average 
temperature, T avg is also given in the table. 

Table 1 Time history of temperature showing how steady state is reached for heat sink with heating coil kept in 

open air for 26 °C ambient temperature. 



Volt- 
age 
(Volts) 


Current 
(Amps) 


Power 
Input 

(W) 


Time 

(s) 


Ti 

(°Q 


T 2 
(°C) 


T 3 
(°C) 


Tavg 

CO 


10 


0.23 


2.3 


0 


26 


26 


26 


26 


65 


28 


28 


28 


28 


289 


30 


30 


31 


30.33 


608 


32 


33 


34 


33 


1034 


33 


34 


35 


34 


1584 


34 


34 


36 


34.66 


2184 


34 


35 


36 


35 


13 


0.31 


4.03 
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32 


32 


32 


32 


122 


33 


34 


34 


33.67 
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35 


35 


35 
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38 


37.33 
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38 
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40 
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41 


41 
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41 


41 
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40 


41 


40.33 


367 


42 


43 


44 


43 


812 


44 


46 


46 


45.33 


1412 


44 


47 


47 
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Fig. 3 shows the plots of the time versus the average temperature, T avg , for different values of the heat 
input. From this figure one can observe the temperature reaching steady state for each case of the heat input. The 
final value of the average temperature T avg , is plotted against the heat input in Fig. 4. 
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Fig. 3 Time history of the average temperature for different values of heat input 
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Fig. 4 Steady state values of average temperature of the heat sink plotted against the corresponding heat input 

values. 

The overall heat transfer co-efficient is computed for the above heat input values from equation (3). 
The values of the heat transfer co-efficient are presented in Table 2 and are plotted against (T avg - T amb ) in Fig. 
5. ' " 

Table 2 Heat transfer coefficient values computed from steady state temperatures obtained from experiments for 
heat sink with heating coil kept in open air (ambient temperature, T a = 26 °C). 



Voltage 


Current 


Heat 


Temperature, 


Heat 


(Volts) 


(Amperes) 


input 


T avg (°C) 


transfer 






(Watts) 




coefficient, 










h 

(W/m 2o C) 


10 


0.23 


2.3 


35 


11.77 


13 


0.31 


4.03 


41.33 


12.11 


16 


0.38 


6.08 


46 


14.00 


19 


044 


8.36 


51.67 


15.00 
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Fig. 5 Heat transfer coefficient values as computed from experimentally obtained steady state temperatures of 

the heat sink. 



IV. THEORETICAL MODELLING 

2-D Analysis 

In 2-D analysis, the 2-D conduction equation with heat generation is solved. The equation is as given 
below [14] 

,d 2 T , d z T. Q _ ... 

Where Q the source term and k is is the thermal conductivity of the material. 

Eq. (4) is used to model the heat transfer in the vertical base plate. The finite volume discretization is 
used to approximate the partial differential equation and the resulting algebraic equations are solved using 
Gauss-Seidal approach. The vertical base plate is divided into a number of finite volume cells. A typical finite 
control volume around a point P of the plate is shown in the Fig. 6. T P represents the temperature at P while T w , 
T E , T s and T N represent the temperatures at the West, East, South and North neighbours of the point P. Ax and 
Ay are the dimensions of the control volume while (8x) w , (5x) e , (8y) s , and (5y) n are the distances of the West, 
East, South and North neighbours from the point P respectively. The finite volume discretization equation can 
be written, following the notation in Patankar [15], as 
a P T P =a E T E + a w T w + a N T N + a s T s + b (5) 
Where 

k e tAy 



k w tAy 
k n tAX 

(«y)n ' 

k s tAy 

as- 

(«y) s 

t = thickness of the plate 



(5a) 
(5b) 
(5c) 
(5d) 




Tig. 6 The control volume 

The source term Q can have the following components 

Q=Qgen+Qflux+Qconv+Qfin (6) 

where 

Qgen =heat generated within the control volume 

= q gen tAx Ay, (6a) 
Qfiux=heat flux through the surface 

= qflux Ax Ay, (6b) 
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Qconv= convec tive heat transfer with the surroundings 

= h face Ax Ay (T P -T amb ), (6c) 
Q fin =heat transferred with the surroundings through fin 

=J?hfm Af ln (T P -T amb ) (6d) 
The coefficient a p and the constant b in Eq. 5 can now be given as, 
a P = a E + a w + a s + a N + h face Ax Ay+ ?7h fm A fin , (7) 

b = q gen tAx Ay+Q flux Ax Ay+ h face Ax Ay T amb + J7h Hn (8) 

In the equations above q gen is the heat generated per unit volume, q flux is the heat flux per unit area, 
h face is the heat transfer coefficient for exchanging heat with surroundings, r\ is fin efficiency, h fin is the 
effective heat transfer coefficient of fin, A fin is the effective surface area of fin through which heat is exchanged 
with surroundings and T amb is the ambient temperature. 

The heat transfer coefficient h tace is modelled using the standard correlations [13], [14]: 
Nu = 0.59 (Gr Pr)° 25 (9) 
Where 

Nusselt number Nu = ^1 (9 a ) 

Grashof number Gr= iM^l (%) 

Prandtl number Pr= (9c) 
Fin efficiency r\ is given by 

tonMml) 
ml 

Where 

mL= L c 15( (-^) 05 (10a) 

L C =L+J (10b) 
A m =t fin L c (10c) 
tf in is thickness of fin. 

The heat transfer coefficient h fin is modelled using the standard correlations [13] [14]: 
for upper surface of fin 

Nu= 0.54 (Gr Pr) 025 (11) 
for lower surface of fin 
Nu= 0.27 (GrPr) 025 (12) 

h fin is obtained from equations (11) and (12) by first obtaining the average Nusselt number as 
Nu= 0.405 (GrPr) 0 25 (13) 

Nu = hfa£s_Le (U) 

The Grashof number in equations (11), (12) and (13) is given by 

Gr=^ ' (14a) 

Pr= ^ (14b) 

where L e is given by L e = — — — (14c) 

2 Qfin + b fin ) 

A fin in equation (6d) is given as 
A fm =2Axl fin (15) 
where 

lf in is length of the fin, 
bf in is breadth of the fin. 
Boundary Conditions: 

Convective heat transfer boundary condition has been applied on all the four edges of the plate viz. 
At the boundary, in general, 

k(^) B =h(T B -T amb ) (16) 

Where T B is the temperature at the boundary, 

(^)b is the temperature gradient at the boundary. 

The resulting algebraic equations are solved using Gauss-Seidal approach as mentioned earlier. Iteration 
method has been used for achieving the convergence of the solution of the equations. A computer programme 
has been written in C language for computation of temperature field. 
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A mesh of size 60x32 has been used for the solution. It is uniform in x-direction and non-uniform in y- 
direction. The details of mesh co-ordinates are given in the appendix. 

V. RESULTS AND DISCUSSION 

The value of heat input and the corresponding heat flux value is given as input to the C-program to get 
the temperature field by using the approach as given in section IV above. Temperature values corresponding to 
the three thermocouple locations are noted from the temperature field and average of these three temperatures, 
T avg is calculated. This is done for different values of the heat input. Heat input values, and the corresponding 
heat flux values, the average temperature of the heat sink and heat transfer coefficient is shown in Table 7. Heat 
transfer coefficient is computed from equation similar to Eq. (3). Fig. 7 shows the average temperature of the 
heat sink plotted against the corresponding heat inputs. Fig. 8 shows the heat transfer coefficient values plotted 
against the corresponding heat inputs. Temperature profile showing the variation of temperature along the mid- 
line (x=60e-3) as function of the vertical distance y is shown in Fig. 9. 

Table 7 Heat sink temperature as evaluated by 2D code for heat sink with heating coil kept in open air 



Heat 


Heat 


Temperature, 


Heat 


input 


flux 


Tavg ( C) 


transfer 


(W) 


(W/m 2 ) 




coefficient, 
h (W/m 2o C) 


2.3 


1288.52 


34.88 


11.93 


4.03 


2257.70 


39.72 


13.54 


6.08 


3406.16 


44.60 


15.06 


8.36 


4683.47 


49.99 


16.06 




Fig. 7 Temperature obtained by 2-D analysis 




Fig. 8 Heat transfer coefficient obtained by 2-D analysis 
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V. COMPARISION OF EXPERIMENTAL RESULTS WITH THEORETICAL RESULTS 

Comparison betwee n heat sink temperatures from experiment and theory is shown in Fig. 10. 
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Fig. 10 Comparison between heat sink temperatures from experiment and theory ■ 
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Comparison between heat transfer coefficient from experiment and theory is shown in Fig. 1 1 
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Fig. 11 Comparison between heat transfer coefficient from experiment and theory ™ h «"p*»''"*n*»i 

VI. SUMMARY AND CONCLUSION 

Temperature of heat sink is measured experimentally. The experimentally measured temperatures have 
been compared with those predicted by the theory and have been found to compare well with each other. 
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APPENDIX 

The mesh used in this work (section IV) of 60x32 i.e. it has 60 divisions in x-direction and 32 divisions in y- 
direction. 

The mesh is uniform in x-direction with Ax — 0.002 m. The y co-ordinates are as shown below: 
y (in m) = 0.000000, 0.001563, 0.003125, 0.004688, 0.006250, 0.007813, 0.009375, 0.010938, 0.012500, 
0.014063, 0.015625, 0.017188, 0.018750, 0.020313, 0.021875, 0.023438, 0.025000, 0.026563, 0.028125, 
0.029688, 0.031250, 0.032813, 0.034375, 0.035937, 0.037500, 0.039062, 0.040625, 0.042187, 0.043750, 
0.045312, 0.046875, 0.048437, 0.050000. 
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